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Abstract. We present the results of a multi-frequency study of the structural evolution of the VLBI jet in the BLLac object 
0716+714 over the last 10 years. We show VLBI images obtained at 5 GHz, 8.4 GHz, 15 GHz and 22 GHz. The milliarcsecond 
source structure is best described by a one-sided core-dominated jet of ~ 10 mas length. Embedded jet components move su- 
perluminally with speeds ranging from 5 c to 16 c (assuming z = 0.3). Such fast superluminal motion is not typical of BLLac 
objects, however it is still in the range of jet speeds typically observed in quasars (10 c to 20 c). In 0716-1-714, younger com- 
ponents that were ejected more recently seem to move systematically slower than the older components. This and a systematic 
position angle variation of the inner (1 mas) portion of the VLBI jet suggests an at least partly geometric origin of the observed 
velocity variations. The observed rapid motion and the derived Lorentz factors are discussed with regard to the rapid Intra-Day 
Variability (IDV) and the y-ray observations, from which very high Doppler factors are inferred. 
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1. Introduction 

■ The S5 blazar 0716+714 is one of the most active BLLac 
[ objects. It is extremely variable on time-scales from hours 
' to months at all observed wavelengths from radio to X-rays. 

The redshift of 0716+714 is not yet known. However, opti- 
cal imaging of t he underlying g alaxy provides a redshift es- 
I timate of z ^ 0.3 dWagner et al.ll9 96). A recent X-ray study by 

■ iKadler et al.l 12004*) suggests a much lower redshift of 0. 1, but 
I since this needs further confirmation a redshift of 0.3 is used 

■ throughout this paper. In th e radio bands 0716+714 is an intra- 
, day v ariable (IDV) source dWitzel et alJ ll986: 'Heesche n et alJ 

119871) . It exhib its a flat radio spec t rum at frequencies up to at 
least 350 GHz. lOuirrenbach et all ( Il99ll) found strongly cor- 
related IDV between radio and optical bands. This and the 
correlated variations between X-ray and optical and the si- 
multaneous variations between optical and radio strongly sug- 
gest an intrinsic origin of the intraday variability dWagner et alJ 
Ll996; Wagner & Witzel 1995). From the observed IDV at cm- 
wavelengths a typic al brightness temperature of Th = lO'^ -^ K 
to lO'^K is derived feaus et alJEooi I Wagner et a l.lll993l) . A 
Doppler factor of 15 to 50 would be needed to bring these 
brightness temperatures down to the inverse-Compton limit of 
lO'^K. 

Very Long Baseline Interferometry (VLBI) studies span- 
ning more than 20 years at cm-wavelengths show a core- 
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dominated evolving jet extending several lOmilliarcseconds 
to the north (Eckart et al. 1986, 1987; Witzel et al. 1988). The 
VLB I Jet is misaligned w ith respect to the VLA jet by ~ 90 ° 
(e.g.. lSaikiaetai]ll987l) . In the literature the jet kinematics 
in 0716+714 are discussed controversially. There exist several 
kinematic scenarios with motions ranging from 0.05 mas yr"' 
to 1 . 1 mas yr"' . However, some of these earlier kinematic mod- 
els are based on only two or three epochs separated by sev- 
eral years, which sometimes could have led to ambiguities 
in the identification of fast components teckartetal. 1987"; 
Witzel et a l. 1988; Schalinski et al. 1992; Gabuzda et al. 1998; 
Perez-Torres et alJl2004l) . More recent studies based on more 
data measured higher jet speeds jTian et al.l200ltljorstad etal] 
2001b; Kellermann et al. 200%, but the results were still incon- 
sistent between the two studies. Thus, it is at present not clear 
if 0716+714 is a superluminal source and and how fast the mo- 
tion of the VLBI jet components is. 

In this paper, we present and discuss our results from a re- 
analysis of the last 10 years of VLBI data on 0716+714, ob- 
tained at frequencies between 5 GHz and 22 GHz. Our analysis 
includes 26 observing epochs listed in Tabled Throughout this 
paper we will use a flat universe, with the following parame- 
ters: A Hubble constant of Ho - 71 kms"' Mpc"', a pressure- 
less matter content of Qm - 0.3 and a cosmological constant of 

- 0.7. With these constants an angular motion of 1 mas/yr 
corresponds to a speed of 18.8 c at z = 0.3 (6.6 c at z = 0.1). 

In Sect. 12] we briefly describe the observations and the data 
reduction, the model fitting process and present the final im- 
ages. In Sect. |3]we present the cross-identification of the indi- 
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vidualjet components and a new kinematic model for the jet in 
0716+714. From this we derive the jet speed and orientation. 
We summarise our results in Sect.|3 

Table 1. Observation log. Listed are the observing epoch, fre- 
quency V, total flux density 5tot(v), beam size, beam position 
angle, peak flux density 5 peak and the lowest contour at 3 cr of 
the map in Figs.^&El 



Epoch 


y 


S tot 


Beam 




^ peak 


3cr 




[GHz] 


[Jy] 


[mas X mas], [°] 


*- heani 


r mJy T 

beam 


1992.73°'" 


5.0 


0.67 


0.8 X 1.0, 





0.61 


1.5 


1992.85" 


22.2 


0.75 


0.9 X 2.7, 


-55 


0.64 


3.0 


1993.71 " 


22.2 


0.40 


0.2 X 0.2, 


35 


0.31 


2.0 


1994.21 


8.4 


0.32 


0.8 X 1.4, 


80 


0.26 


0.7 


1994.21 


22.2 


0.34 


0.3 X 0.5, 


79 


0.30 


1.8 


1994.67 


15.3 


0.46 


0.5 X 0.7, 


-29 


0.40 


1.8 


1994.70 


5.0 


0.36 


0.8 X 1.2, 


-15 


0.29 


2.7 


1995.15*'' 


22.2 


0.75 


0.3 X 0.5, 


-3 


0.71 


4.0 


1995.31 


22.2 


0.43 


0.3 X 0.4, 


39 


0.40 


4.2 


1995.47 


22.2 


0.20 




7 


0.18 


3.0 


1995.65 


8.4 


0.31 


0.9 X 1.0, 


54 


0.26 


0.4 


1995.65 


22.2 


0.34 


0.3 X 0.4, 


-40 


0.28 


1.8 


1996.34 <3) 


22.2 


0.27 


0.3 X 0.9, 


48 


0.22 


3.2 


1996.53 


15.3 


0.26 


0.5 X 0.7, 


-15 


0.22 


0.7 


1996.60 <3) 


22.2 


0.25 


0.4 X 0.4, 


9 


0.21 


3.0 


1996.63 


5.0 


0.22 


1.4 X 2.2, 


41 


0.18 


1.3 


1996.82 <2) 


15.3 


0.26 


0.4 X 0.7, 


21 


0.25 


2.0 


1996.90 <3> 


22.2 


0.29 


0.4 X 0.5, 


-35 


0.27 


2.0 


1997.03 


8.4 


0.21 


1.1 X 1.4, 


7 


0.18 


0.0 


1997.58 <3) 


22.2 


0.98 


0.4 X 0.5, 


-33 


0.90 


4.0 


1997.93 


8.4 


0.43 


0.6 X 1.6, 


33 


0.38 


0.5 


1999.41 


8.4 


1.02 


0.5 X 1.5, 


-7 


0.92 


0.5 


1999.55 <2' 


15.3 


1.25 


0.5 X 0.8, 


4 


1.14 


1.2 


1999.89 "•<'> 


5.0 


0.63 


1.8x2.4, 


-60 


0.54 


0.7 


2000.82* 


5.0 


0.55 


0.8 X 0.9, 


-28 


0.48 


0.1 


2001.17 <2) 


15.3 


0.65 


0.4x0.8, 


31 


0.54 


0.9 



Note: The array used was the VLBA, unless indicated by a footnote. 
Epochs in bold face denote ow n data, a: Global array, b: VLBA+Eb 
Referenc es: f nilBritzen et aljfl 999.1 & in prep . (20 05), (2): 
FKellerm ann et alj <1998h &IZensus et S" J2002h. f3):IJorstad et alj 
l2001bl) .74i lFev & CharlotH2000l) . ('51: Iros et alji200ll) 
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2. Observations and data analyses 

All data were correlated in the standard manner using the 
MK III correlator in Bonn and the VLBA correlator in Socorro. 
Part of the post-correlation analysis was done using NRAO's 
Astronomical Image Proces sing System (AiPs) and the C altech 
VLBI Analysis Programs dPearson & Readheadlll984l) . Data 
from other observers were provided as amplitude-calibrated 
and fringe fitted (u, v) FITS -files. The imaging of the source 
including phase and am plitude self-calib ration was done using 
the CLEAN algorithm ('Hosbom"l974") and SELFCAL pro- 
cedures in DiFMAP ( Shepherd 1994 ). The self-calibration was 
done in steps of several phase-caUbrations followed by careful 
amplitude calibration. During the iteration process the solution 
interval of the amplitude self-calibration was shortened from 



intervals as long as the whole observational time down to min- 
utes. The resulting maps are presented in Figs.n&|2] Here, the 
jet is clearly visible to the north at PA. ~ 15 °, and is slightly 
bent. At the lower frequencies we can follow the jet up to a dis- 
tance of 10 mas to 15 mas from the core, whereas at the higher 
frequencies the jet is visible up to 3 mas. 

2.1. Model fits 

After imaging in Difmap we fitted circular Gaussian compo- 
nents to the self-calibrated data to parameterise the source 
structure seen in the VLBI maps. We used only circular com- 
ponents to reduce the number of free parameters and to sim- 
plify the analysis. A summary of the Gaussian component pa- 
rameters from this model fitting is given in Table |2l For all 
model fits we used the brightest component as a reference and 
fixed its position to (0,0). The positions of the other compo- 
nents were measured relative to this component, which we as- 
sumed to be stationary (see Sect. l3.1l for details). A maximum 
uncertainty of 15 % in the flux density was estimated from the 
uncertainties of the amplitude calibration and from the formal 
errors of the model fits. The position error is given by Ar - ^ 
(Fomalont 1989), where cr is the residual noise of the map after 
the subtraction of the model, the full width at half maximum 
(FWHM) of the component and 5 peak the peak flux density. 
This formula tends to underestimate the error, if the peak flux 
density is very high or the width of the component is small. 
Therefore we included an additional error arising from the posi- 
tion variations during the model fitting procedure. In Fig.|3 the 
core distances of individual VLBI components derived from 
our model fits at various frequencies are plotted against time. 
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Fig. 3. Core separation as a function of time for the individ- 
ual modelfit components. Data from all frequencies are com- 
bined. Possible frequency-dependent position shifts are less 
than 0. 1 mas and are not corrected. The solid lines show the Un- 
ear fits to the path for each component. To show more clearly 
the well-defined components we do not show two of the far- 
thest data points of CI at r = 11.2 mas and r - 11.9 mas and 
one data point of C2 at r = 6.9 mas. 
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Fig. 1. All contour maps of 0716+714 at 5 GHz (top row), 8 GHz (middle) and 15 GHz (bottom). The maps are convolved with 
circular beams of 1 .2 mas at 5 GHz, 0.8 mas at 8 GHz and 0.5 mas at 15 GHz. Total flux density, original beam size and the level 
of the lowest contour at three cr are given in the observing log (Table^l- Contours are increasing by steps of two. 



3. Results and discussion 

3.1. Identification of thie components 

To investigate the kinematics in the jet of 0716+714, we cross- 
identified individual model components along the jet using 
their distance from the VLBl core, flux density and size. Since 
the observations were not phase-referenced, the absolute po- 
sition information was lost, and it is impossible to tell which 
components, if any, are stationary. To register the images and 
to test the stability of the VLBI-core position, we compared the 
core separations of the jet components of each epoch with re- 
spect to the adjacent observations, but we could not find any 
systematic position off'sets. We also checked for position shifts 
between the jet components at higher and lower frequencies 
due to opacity effects (e.g., L obanov 1998). Inspection of close 
or simultaneously observed frequency pairs does not show a 
systematic frequency shift between higher and lower frequen- 
cies larger than 0.1 mas (position accuracy between 15 GHz 
and 22 GHz is < 0.1 mas). Therefore, we conclude that the core 
position is stationary and that all component positions can be 
measured with respect to this brightest component. 



The cross identification of moving VLBI components be- 
tween different times and frequencies is not unambiguous and 
depends on the dynamic range of the individual maps and on 
the time sampling. When we started our analysis we used the 
following scenarios as working hypotheses: 

1 . Stationary core with stationary or slow-moving jet, 

2. Stationary core with fast moving jet, 

3. Non-stationary core with fast or slow-moving jet, 

4. Non-monotonic motion of the jet. 

As more data became available, we could rule out most of 
these identification schemes, and we were left with a scheme 
that assumes a stationary core and relatively fast component 
motion. The different schemes were tested, first separately at 
each frequency and later with all modelfits combined. 

Supported by a graphical analysis, which is presented in 
Fig. 13 we could obtain a satisfactory scenario for the kinemat- 
ics in the jet of 0716+714. This scenario consists of 1 1 identi- 
fied superluminal components moving away from the core. As 
an example of the procedure used for all components we show 
the identification scheme for the separate frequencies for com- 
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Fig. 2. All contour maps of 0716+714 at 22 GHz. The maps are convolved with a circular beam of 0.3 mas. Total flux density, 
original beam size and the level of the lowest contour at three cr are given in the observing log (Table^. Contours are increasing 
by steps of two. 



ponent C5, CI, C8 and C9 in Fig.|3 Again the frequency shifts 
between the trajectories are typically smaller than 0. 1 mas to 
0.2 mas and are smaller than the measurement uncertainty Ar. 

Table |3l summarises the angular separation rate fi of the 
VLBI components derived from the linear fits of r(f) and the 
corresponding apparent speed jSapp - Vapp/c for an assumed 
redshiftof 0.3. 

3.2. The kinematic model 

The components in this scenario move with 0.2 mas yr"' to 
0.6 mas yr ' in the inner part of the VLBI jet (r < 3 mas) and 
with up to 0.9 mas yr"' in the outer regions. Due to the inho- 



mogeneous time sampling, it is sometimes difficult to iden- 
tify components over a long time interval. However, especially 
between 1993 and 1996, when we made many measurements 
spaced by only a few months, our proposed component identi- 
fication gives the simplest and most reasonable fit to the data. 
Since most of our data were obtained at high frequencies, 
where the outer region of the jet is faint and already partly re- 
solved by the interferometer beam, the parameter of the corre- 
sponding jet components (at r > 4 mas) are not as well con- 
strained as the inner jet components and therefore tiiey have 
larger positional errors. Despite this, there is still a very dis- 
tinct tendency that the older components (CI - C3), which now 
are located at large core separations, move faster than the com- 
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Fig. 4. Core separation as a function of time for the components 
C5, C7, C8 and C9 separated for the different frequencies. The 
soHd Unes show the Unear fits to the path for each component. 
Shifts between the paths at dififerent frequencies are visible, but 
are typically smaller than the measured uncertainties. 

ponents located in the inner jet (C4 - CIO). This behaviour is 
also visible in Table|3 which shows a clear trend with system- 
atically decreasing speeds between component CI (16 c) and 
component CIO (4.5 c). The coarse time sampling for the in- 
dividual jet components unfortunately does not allow us to fit 
acceleration to r(t). Future and more densely sampled VLBI 
observations are required to show whether the inner jet compo- 
nents indeed move linearly. 

Table 3. Proper motions in the jet of 0716-1-714. The number of 
data points which were used for the fit for each component is 
given in column 2. The last column gives the back-extrapolated 
ejection dates of the components, which result from the linear 
fits. 



Id. 


# 


/J [mas/yr] 




= 0.3) 


Ejection date 


CI 


4 


0.86 ±0.13 


16.13 + 


2.36 


1986.7 


± 1.3 


C2 


5 


0.63 + 0.10 


11.87 + 


1.82 


1988.5 


± 1.0 


C3 


8 


0.62 ± 0.04 


11.59 + 


0.76 


1989.7 


±0.3 


C4 


7 


0.44 ± 0.01 


8.27 + 


0.24 


1989.8 


±0.1 


C5 


7 


0.43 ± 0.02 


8.16 + 


0.29 


1990.8 


±0.2 


C6 


9 


0.41 ± 0.02 


7.80 ± 


0.33 


1992.3 


±0.1 


C7 


17 


0.37 + 0.01 


6.89 ± 


0.19 


1993.2 


±0.1 


C8 


11 


0.32 ± 0.01 


6.04 ± 


0.20 


1994.4 


±0.1 


C9 


10 


0.26 ± 0.01 


4.98 ± 


0.27 


1995.1 


±0.1 


CIO 


3 


0.24 + 0.01 


4.52 ± 


0.45 


1996.5 


±0.2 


Cll 


2 


0.29 + 0.02 


5.48 ± 


0.55 


1998.2 


±0.3 



The luminosity distance dh was calculated adopting the follow- 
ing relation 



77(1, Qm) - ri( 



1 



1 +z 



n„) 



(1) 



where 

ri{a,Q.^) = 2 Vs3 -H 1 



1 s 

^ -0.1540-7 -I- 0.4304-5- 



H-0. 19097- 



1 -Qn 



and 



■ 0.06694 L?" 
1 



with 



which is an analytical fit (Pen 1 9991) to the luminosity dis- 
tance for flat cosmologies with a cosmological constant. For 
the range of 0.2 < Qm < 1 0, the relative error of the distance 
is less than 0.4 % for any given redshift. Using 



^app — 



C(l -^z) 



(2) 



for the transformation of apparen t angular separation rates 
(//) into spatial apparent speeds (/Sapp) rP earson & ZensusI 19871) 
and assuming a redshift of 0.3 for 0716H-714, one milliarc- 
second corresponds to 4.4 pc. The measured angular separa- 
tion rates then correspond to speeds of 4.5 c to 16.1 c. For the 
components C3 to C9, for which the proper motion determi- 
nation is the most confident (over seven data points), the aver- 
age motion is /i = (0.41 + 0.11)masyr"' which corresponds to 
ySapp = 7.7 + 2.1. A preliminary analysis by Tian et al. (200ll). 
who used only part of the data presented here, gave very simi- 
lar speeds. A recently published kinematic analysis of the data 
from the VLBA 2 cm Survey a lso measured motion o f 10 c 
to 12 c in the jet of 0716+714 jKeUermann et alJl2004l) . The 
speeds found in our paper are lower than the 16 c to 21 c found 
by Jorstad et al. ( 2001b) '. who observed the source over only 
a short time interval (3 yr). The difference between their and 
our results may be explained by slightly different (and unfor- 
tunately not unambiguous) map parameterisations using differ- 
ent Gaussian models or by non-linear component motion, with 
possible acceleration for shorter periods of time. 

With an average apparent jet speed of about 8 c and 
likely higher speeds of up to 15 c to 20 c, 0716+714 dis- 
plays considerably faster motion than other BLLac objects, 
for w hich speeds of < 5 c are regarded as no r mal (e.g., 
iVermeulen & CoherJ 1 1994 lOabuzda et all i200d IRos et alJ 
I2OO2I) . In this context it appears that 0716+714 is an extreme 
BLLac object with a jet speed (and Lorentz factor) much 
higher than regarded typical for BLLac objects and close to 
the speeds of 10 c to 20 c seen in quasars. 

3.3. Kinematics and geometry of the jet 

Using the measured motion, (11.6 + 0.8) c, of C3, the fastest, 
best-constrained jet component in our model, we can place lim- 
its on the jet speed and orientation of 0716+174. Adopting 



app 



/3smd 
1 -yScose 



(3) 



we find that the jet inclination that maximises the apparent 
speed, flmax is 4.9 °. From that, one derives a minimum Lorentz 
factor, 7niin = (1 + ySapp)'^^ 11.6, which corresponds to a 



' Jorstad et al. originally published 11c to 15 c using Hq = 
100 kms"' Mpc"' and qo = 0.1. We corrected these numbers for the 
cosmological parameters used in this paper. 
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Doppler factor of 6, 



„(l-AninCos6lmax) ' = 11.7,where 



1 - 7 . For smaller viewing angles (6 



0) the 

Doppler factor approaches 6 - 2y, which yields 6 = 11 .7 to 
23.4 at viewing angles between 4.9° a nd 0° for componen t 
C3. Including the higher values found bv*Jorst ad et al 
and the highest speed derived from our data (16.1 c), a Doppler 
boosting factor of up to 5 > 40 appears possible. 
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Fig. 5. The Doppler factor versus the apparent speed for con- 
stant intrinsic Lorentz factor, y, (solid lines) and constant view- 
ing angles Q (dotted lines). The range of measured speeds is in- 
dicated by the dashed lines and the grey shaded area marks the 
possible values for j and Q (see text for details). 



At cm-wavelengths, 0716+714 is a prominent intraday 
variable source, which shows amplitude variations of 5 % 
to 20 % on time-s cales of 0.25 d to 2d ( Wagner et al. 1996: 
iKraus et alJ l2003l). Correlated radio-optical IDV observed 
in this source ('Ouirrenbac h etalJll99l[ IWagner etalJ ll996^. 
[oian et al. 1996) suggests that at least some fraction of the ob- 
served rapid variability has an intrinsic source origin and can- 
not be attributed to refractive interstellar scintillation (RISS) 
alone , as is done for some other IDV sources (e.g.. ,Rickett 
l200lt lOian et al.ll200lt Ikedziora-Chudczer et alJl200lh The 
presence of the observed broad-band correlations of the vari- 
ability and recently detected IDV at 9 mm wavelength, where 
RIS S should not play a dominant role due to its X"'^ dependence 
(cf . iKrichbaum etal]l2002t iKraus et alj bool . further support 
a non negligible intrinsic contribution to the IDV in 07 16+7 14. 
In the following we use an average typical brightness temper- 
ature derived from the IDV observed at 6 cm of 10'^^ K to 
10'^ K. Slightly higher values of up to a few times 10"^ K were 
measured only occasionally. To bring this brightness temper- 
atures down to the inverse Compton limit of 10'^ K, Doppler 
factors in the range of 5 = (Tb/IO'^K)" a; 15 to 50 are re- 
quired. Adopting these Doppler factors, we obtain jet Lorentz 
factors of 7 a; 8 to 25 and viewing angles of < 1 °. 

The relation between the Doppler factor, intrinsic Lorentz 
factor, viewing angle and apparent speed is illustrated in Fig.|5] 
The Doppler factors derived from IDV lie within the grey 
shaded area. 



If we now try to explain the observed range of apparent 
component speeds through variations of the viewing angle (un- 
der the assumption of a constant Lorentz factor along the jet), 
we cannot reach slow apparent speeds of ~ 4.5 c without vi- 
olating the lower limit of the Doppler factor of 2.1 from syn - 
chrotron self-Compton (SSC) models (Ghi sellini et alJll993l) . 
Therefore, we can exclude Q > 5° and the only solution re- 
maining is to decrease the viewing angle which automatically 
leads to higher Doppler factors. The grey shaded area in Fig.|5] 
marks the allowed region of Lorentz factors and viewing an- 
gles. It becomes obvious that we need at least a Lorentz factor 
of 7 a 15 and a viewing angle of the VLBI jet of 6* ~ 2 ° to ex- 
plain the large range of observed apparent speeds as an effect 
of spatial jet bending. These values are consistent with those 
derived from IDV. 

3.4. F/ux density evolution 

In Fig.|6lwe show the evolution of the flux density of VLBI 
components at 22 GHz with their increasing separation from 
the core. We show this frequency because a large fraction of the 
data were obtained at 22 GHz (11 epochs). With the exception 
of a few data points, the components fade as they travel down 
the jet. This agrees qual itatively with the theory of a conically 
expanding jet (e.g., Blandford & K6nigllll979l) . As the VLBI 
components expand they become optically thin and their spec- 
tra steepen. This is illustrated for component C7 in Fig.0 The 
spectral index a is defined as 5 oc v". Knowing the turnover 
frequency and the size of a jet component allows u s to calcu- 
late th e strength of the magnetic field in the jet (e.g.. lMarscheJ 
Em as 

B = 10-5 b(a) f''^ gauss, (4) 

Smi^ +Z) 

where b(a) is a t abulated parame ter dependent on the spectral 
index (Tab. 1 in iMarsch er' 1983) with a value ranging from 
1.8 to 3.8 for optically thin emission (a = -0.25 to -1.0), 
Sm (in Jy) is the flux density at the turnover frequency Vm 
(in GHz) and 6 (in mas) is the size of the component. If we 
calculate B for each of our modelfit components, without cor- 
recting for the Doppler factor 6, we get the apparent magnetic 
field along the jet. Since we do not have adequate spectra for 
all positions along the jet we will consider a decrease of the 
turnover frequency along the jet proportional to H', to correct 
for the expansion of the component. The best results are ob- 
tained when we start with a turnover frequency of 4 GHz for the 
inner jet (r a 0.2 mas), that is plausible from the early spectrum 
(1994) of component C7 (Fig.0 though the turnover is not ac- 
tually observed, and use p - -1.3, which is in good agree - 
ment with the relativistic jet model ( lMarscher&Geailll985l) . 
For all components along the jet we derive magnetic fields in 
the range of lO ^^^G to 10"^ 5 G. Since the minimum mag- 
netic field strength that should be present due the energy den- 
sity o f the cosmic m icrowave background is a few microGauss 
(e.g.. ivan der Laan & Pero la 1969) the 10"^ 6 G already implies 
a lower limit of the Doppler factor of about 10. At this stage the 
core component is excluded from the analysis because both the 
source size and the turnover frequency are poorly known. 
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Fig. 6. Flux density of the VLBI components at 22 GHz plot- 
ted versus core separation. The components fade as they travel 
down the jet. 



Assuming equipartition between the energy of the particles, 
Eg, and the energy of the magnetic field, Eb, one can also derive 
the minimum magnetic field from the synchrotron luminosity, 
L. 



L ^ 4ndl(l+z) f 
Ee = f(a,vuV2) L B 



S dv 



1.5 



(5) 
(6) 



where f(a, vi , V2) is a tabulated function (e.g. IPacholczv^ 
Il97nl) and vi and V2 are the upper and lower cutoff frequencies 
of the synchrotron spectrum with typically values of 10^ Hz to 
10" Hz. Taking the energy density of the magnetic field to be 
I/SttB^ and assuming spherical symmetry the total energy of 
the source is 



(l+k)E, + Es 

(1 +k)f(a,vuV2)LB-^-^ + \r^B^-, 

6 



(7) 



where k is the energy ratio between the heavy particles and the 
electrons and R is the size of the component. The ratio k de- 
pends on the mechanism of generation of the relativistic elec- 
trons, which is unknown at the present. It can range from k ^ \ 
for an electron-positron plasma up to 2000, if most of the 
energy is carried in the protons. We used k « 100 i n our calcu- 
latio ns, which seems to be a reasonable value (e.g.. |Pacholczvkl 
Given that E^cc B ' ^ and E^ B^, Etot has a minimum 
and we find 

/ 9 \^'^ 

= 5.37 -10'^ (SmVmdlR-Y' ■ 

with k = 100, /(-0.5, 10^ 10") = 1.6 ■ 10^ and 5„ in Jy, 
Vni in GHz and R in cm. Using z = 0.3 and the same values 
for Sm, Vm and 6 as in Equation |4] we obtain magnetic fields 
between 10 "* G and 10"^ G for the jet. These values are larger 
than those derived before, which means that either the particle 
energy dominates and we do not have equipartition or that the 



emission is relativistically beamed. In the latter case, the dif- 
ferent dependence of the magnetic field on the Doppler factor 
in Equation|3(Z? oc S) and in Equation|8](Bmin ^ ^^/Vck+i-j ^^jj 
be used to calculate the Doppler factor from B^i„/B = d^^''". 
Comparing the results from all jet components yields Doppler 
factors of 4 to 15 with a mean value of 9. These values are 
slightly lower than those derived from the kinematics, but given 
that the calculations are based on many assumptions the values 
agree reasonably. 
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Fig. 7. Spectral evolution of component C7. 

In Fig.|S]we show the long-term radio variability and the 
spectral index of 0716+714 over the time range in which the 
VLBI components were born. We show a combined data set 
with flux density measurements at 5 GHz and 15 GHz from 
the U MRAO flux density monitoring program dAller et alJ 
I2OO3I) and from our flux-density monitoring performed with 
the 100 m radio telescope at Eff'elsberg (Peng et aL .2000.) . A 
detailed discussion of the flux den sity variability using these 
and other data was recently done bv Ra iteri et alJ ll2003l) Here, 
we restrict the discussion to a possible correlation of the radio 
variability with the ejection of new jet components. The ejec- 
tion dates of new VLBI components are indicated by the grey 
shadowed areas. Their widths correspond to the uncertainty in 
the back-extrapolated ejection dates (see Table Although 
these uncertainties are rather large and the time sampling of 
the light curves is not always dense enough, a weak correlation 
between the ejection of new components and the flares in the 
radio bands is obvious. Each of the shaded areas either lies at 
or shortly before the time of a flux density increase of at least 
one of the two observing bands. 

The two outbursts in mid 1992 and early 1995, which are 
best seen in the 2 cm band, are both surrounded by two new 
components, but the time sampling is too sparse to draw any 
further conclusions. In the lower panel of Fig.|S]we calculated 
the spectral index between 5 GHz and 15 GHz from the two 
light curves. Since the light curves were not measured simul- 
taneously we searched for the closest pairs in the data set and 
calculated the spectral index from these pairs. The mean sepa- 
ration between two points is about eight days. The graph itself 
represents an eight-point running average to reduce the noise in 
the data. This figure shows an obvious correlation between the 
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1988 1990 1992 1994 1996 1998 2000 
epoch 

Fig. 8. The long term flux density variability of 0716+714 as 
measured at Effelsberg and with the Michigan radio telescope 
at 5 GHz (top) and 15 GHz (middle) and the spectral index (bot- 
tom), as/isGHz- The shaded areas mark the ejection dates of the 
VLBI components and their uncertainties. 



spectral index variations and the ejection of a new VLBI com- 
ponent. Seven out of eight ejections (except C9) of new VLBI 
components are accompanied by a flattening of the source spec- 
trum. The flattening is in good agreement with an expanding 
component that becomes optically thin first at the higher fre- 
quenc y and later at the lower frequency (^g.- lMarscher & GeaJ 
Il985l) . 

We note that at the times of particularly dense time sam- 
pling (after 1994) a number of minor radio flares are visible, 
which are not related to the ejection of any of the known VLBI 
components. According to the light-house model and other 
related helical jet models (e.g., Camenzind & KrockenbergeJ 
[l992: Roland et al. 1994) it is possible that initial outbursts, 
which are related to the ejection of new jet components are 
followed by secondary flux density variations, which are not or 
only are indirectly related to the component ejection. Of course, 
it is also possible that we missed some jet components in our 
infrequently sampled VLBI monitoring. 

To test the significance of the correlations found by eye, we 
placed eight fields with a similar width randomly over the light 
curves between 1989 and 1999 and determined how many are 
coincident with a flare in the light curve or a flattening of the 



spectrum. For the light curves this test yielded a 50 % prob- 
ability of our results occurring at random. The probability of 
measuring a flattening of the spectrum in seven out of eight 
components selected at random from our spectrum is 5.7%. 
Generally a result is excepted to be statically significant if the 
probability is below 5 %, which means that the correlations of 
the light curve with the ejection dates are likely to be a coin- 
cidence, but given that we might have missed a component the 
correlation of the spectral index is more significant. A coordi- 
nated and much denser sampled multi-frequency flux density 
and VLBI monitoring is necessary to study in more detail such 
outburst-ejection relations. 

3.5. Connection with y-ray flares 

0716+714 was detected at y-rays by the EGRET de- 
tector on board the Compton Gamm a Ray Observatory 



dHartman et al]ll999l:lMattox et al.ll 19971) . Since it is proposed 
that the y-ray emission, like the radio emission, is Doppler 
boosted (iDermer & Schlickeiser 1994) . the y-ray detection of 
0716+714 yields further evidence for a high Lorentz factor in 
the jet. 

Recent studies about the connection between y-ray sources 
and superluminal VLBI-components suggest that these sources 
have higher aver age jet sp eeds than sources which are not de - 
tected at y -ravs jjo rstad et al. |,2001a: Kellermann et alJl2004l) . 
Moreover, Ijorstad e7am2001al) found a correlation between 
the ejection of new VLBI components and the appearance of 
y-ray flares, suggesting that the radio and the y-ray emission 
originate within the same shocked area in the relativistic jet. 
The authors report a y-ray flare of 0716+714 in 1992.2 which 
appears to be coincident with a peak in the linear polarization 
light-curve at 15 GHz. Unfortunately, their 0716+714 VLBI 
data does not cover this period, but with our larger sample 
we found a new VLBI component (C6) with an extrapolated 
ejection date of 1992.3 ±0.1 which is in good agreement with 
the reported y-ray flare. Although we register three additional 
ejections up to 1996 and 0716+714 was observed several times 
by EGRET in this period, no more flares were detected. But 
given that the y-ray measurements have large uncertainties and 
only 18 data points were recorded during these four years, 
this seems not exceptional. Simil ar outburst ejection ev ents 
were also foun d in e.g, 05 28+134 ( Krichbaum et al."l995') and 
0836+710 COtterbein et al. 1998) . Assuming that the correla- 
tion between the new component, C6, and the 1992.2 flare is 
real, this supports the idea that the y-ray emission originates 
from the same region from w hich we observe the radio emis- 
sion re.g..lSikoraet"^lll994l)V 



3.6. Variation of tine ejection position angle and core 
flux density 

The visual inspection of position angle of the innermost portion 
of the VLBI jet in the maps of Figs.^&IJlindicates a variation 
of the R A. of the jet near the core with time. In order to quantify 
this, we used the Gaussian models and fitted a straight line to 
the first milliarcsecond of the jet. For each observing epoch. 
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the position angle of this line provides a good estimate of the 
direction of the inner portion of the VLBI jet. In Fig. |9] (top 
panel) we plot the P.A. grouped in one year time bins. 
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Fig. 9. Top: Position angle of the central 1 mas of the VLBI jet 
grouped in one year time bins. The dashed line is a sinusoidal fit 
to the data. The position angle of the inner jet (ejection angle) 
varies with a period of ~ 7.4 years. Bottom: One year bins 
of the 22 GHz core flux density. The dashed line represents a 
sinusoidal fit with a fixed period of 7.4 years. 



A sinusoidal fit to the data reveals a period of (7.4 + 1 .5) yr. 
The fit yields a reduced of 1.55, which is much better than 
flie Xi of 4.28, obtained for a Unear fit. Unfortunately our data 
cover only one period of the variability and it will be interest- 
ing to see if the trend holds in the future. A possible explanation 
for such periodic variation would be a precessing footpoint of 
the jet. This should also result in a flux density variation of 
the VLBI core, and indeed, there is a tendency that the core 
flux density at 22 GHz is consistent with a periodic variation of 
~ 7 yr (see Fig.|5] bottom panel). Recentlv Rai teri et alJ ( l2003l) 
published an analysis of the optical and radio flux density vari- 
abiUty of 0716+714 and found a periodicity of 5.5 yr to 6.0yr 
for the flux-density variations in the radio regime. This is com- 
patible within the error bars with the periodicity of the P.A. in 
the inner region of the jet. Based on the geometrical consider- 
ations in Sect. 13.31 the apparent peak to peak variation of ss 7 ° 
of the ejection angle corresponds to a smaU change of only 
7 ° sin e a; 0.6 ° of the jet direction. 

Such small variations of a few degrees in the rest frame 
of the source in several years have recently been detected 
also in other AGN, e.g. 3C 34 5 (.Biretta et al. .1986). 3 C 279 
llCarrara et alJ Il993l). 3C 27 3 llKrichbaum etal.l boOll) and 
BLLac ("Stir ling et allEooi ). The most common explanation 
is a precessing jet, which c ould be caused by a binary super- 
massi ve black hole system dBiretta e t alj|l986t iHummel et alJ 
I1992I) or a warped accretion disc ('Lai"2003' and ref. therein). 
However, hydr odynamic eff'ects al so can lead to a precession 
of the jet (e.g.. lHughes et al]l2002l) . In this case the precession 
is introduced by a potentially strong oblique internal shock that 
arises from asymmetric perturbation of the flow. 



3.7. Redshift dependence of the results 

Since all the derived parameters of the kinematics and geome- 
try scale with the distance to 0716+714, knowing the redshift 
is of great importance. The distance to 0716+714 and there- 
with the speeds of the jet components scale nearly linearly 
with z at these low redshifts. Thus a redshift decrease by a fac- 
tor of three would also decrease the Doppler factor by three. 
The dependence of brightness temperature derived from IDV 

is T\, oc (-[^) (e.g., a decrease from z = 0.3 to z = 0.1 would 
reduce T\, by 6.4). Since, oc this would increase the dis- 
crepancy between the Doppler factor, which is needed to re- 
duce Tb and the Doppler factor derived from the kinematics. An 
increase in redshift would bring the derived parameters closer 
together 

Our assumption of the redshift of z > 0.3 is based on the 
non-de tection of an underlying host galaxy by i Wagner et alJ 
('1996). BLLac host galaxies have typical absolute magni- 
tudes of Mr = a -24 and an effective radius of ~ lOkpc 
(Jan nuzi et alJ 1 19971 and ref. therein). The lower limit of 
Wagner et al. corresponds to an absolute magnitude of Mr - 
-20.4, which would be an unusually faint BLLac host and so 
is a very conservative lower limit. An upper limit for the red- 
shift can be given from the lack of absorption by foreground 
galaxies, which results in a redshift of z < 0.5 (Wagner priv. 
comm.). At present it is unclea r if the tentative X- ray detection 
of a spectral line near 5.8keV( lKadier et al. 2004) can be iden- 
tified with the FcKq, line. If so, the redshift of the line would 
indicate a distance of z - 0.10 ± 0.04, much closer than ex- 
pected from the optical measurements. Possible ways out of 
this difficulty could invoke a blue shift of the line or some ex- 
ternal emission process not related to the nucleus of 0716+71. 
A solid confirmation of the X-ray line in future observations is 
urgently needed. 

4. Conclusions 

We analysed 26 epochs of VLBI data at 4.9 GHz, 8.4 GHz, 
15.3 GHz, and 22.2 GHz observed over 10 yr, between 1992 
and 2001, and derived a new kinematic scenario for the jet in 
0716+714. In this scenario the components move with appar- 
ent speeds of 4.5 c to 16.1 c (z - 0.3). These speeds are atyp- 
ically fast for BLLac objects, which typically exhibit speeds 
of < 5 c (e.g..lVermeulen & Cohenll994[ iGabuzda et"aDl200(i 
iRos et alJ ^OOzT Since this analysis is based on much denser 
sampled VLBI data than previous studies, we are convinced 
that we can rule out the so mewhat slower scenarios which were 
presented previously (e.g., W itzel et a i1'l988''Schali nski etalJ 
J992; GabuzdiLeLal, 1998 ) . We find a periodic variation of the 
P.A. of the innermost jet which could be a sign of jet preces- 
sion on a time-scale of (7.4 + 1.5) yr. If this can be confirmed 
in futur e observations the higher speeds found bv .lorstad et alJ 
(l2001bl) might be explained by non-linear motion of the jet 
components. 

No correlation was found between the component ejec- 
tion and radio flux density flares at the cm-wavelengths. There 
seems to be a weak correlation between the ejection of new 
components and the flattening of the radio spectrum. For one 
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of the new components (C6) we could find a reported y-ray 
flare and if the ejection and the flare are really connected this 
would support the idea that the y-ray emission originates from 
the same region from which we observe the radio emission. 

From the component motion in the jet, we obtain a lower 
limit for the Lorentz factor of 11.6 and a maximum angle to 
the hne of sight of 4.9 °. To explain the large range of observed 
apparent speeds as an efl'ect of spatial jet bending, a Lorentz 
factor of y > 15 and a viewing angle of the VLBI jet of < 2 ° 
are more hkely. Under these circumstances, the Doppler fac- 
tor would be (5 « 20 to 30. Such high Doppler factors are 
indeed required to explain the high apparent brightness tem- 
peratures of up to lO'^ K derived from intraday variability at 
cm-wavelengths. 
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Table 2. Results from Gaussian Model fitting and component 
parameters. Speak is the peak flux density, r and are the dis- 
tance and the PA. measured from the core and Q is the FWHM 
of the Gaussian component. 



Epoch 


V 


^ peak 


r 




e 


Id." 




[GHz] 




[mas] 


[°] 


[mas] 




1992.73 


5.0 


596 ± 89.4 


_ 


_ 


0.10 ±0.02 


Core 






37 ±5.5 


0.83±0.12 


17 ±8 


0.14 ±0.03 


C5 






19 ±2.8 


1.92±0.13 


8±3 


0.60 ±0.1 2 


C3 






8 ±1.2 


5.02 ±0.24 


I7±2 


0.93 + 0.19 


CI 






7 ±1.0 


11.46 ±0.49 


13±2 


3.35 ±0.67 


X 


1992.85 


22.2 


708 ± 106.2 


_ 


_ 


0.18 ±0.04 


Core 






12±1.8 


1.37 ±0.06 


12±2 


0.10 ±0.02 


C4 






16 ±2.4 


2.86 ±0.22 


0±4 


1.48 ±0.30 


C2 


1993.71 


22.2 


272 ±40.8 






0.05 ±0.01 


Core 






79± 11.8 


0.19 ±0.02 


15±5 


0.05 ±0.01 


C7 






14±2.1 


0.58 ±0.07 


14 ±6 


0.12 ±0.02 


C6 






14±2.1 


1.22 ±0.09 


20±4 


0.23 ±0.05 


C5 




22.2 


309 ±46.3 


_ 


_ 


0.09 + 0.02 


Core 






15 ±1.9 


0.28 ±0.09 


28 ±18 


0.01+0.00 


C7 






13 ±1.9 


0.80±0.13 


15±9 


0.42 ±0.08 


C6 






7 ±1.0 


1.79±0.16 


14±5 


0.38 ±0.08 


C4 


1994.21 


8.4 


266 ±39.9 


_ 




0.14 ±0.03 


Core 






25 ±3.7 


0.82 + 0.09 


13±6 


0.42 ±0.08 


C6 






10± 1.5 


1.99±0.19 


17±5 


0.76±0.15 


C4 






9± 1.3 


3.27 ±0.22 


12±3 


0.84±0.17 


C2 






4 ±0.6 


4.70 ±0.43 


14±5 


1.48 ±0.30 


X 


1994.67 


15.3 


413±61.9 


_ 


_ 


0.07 ±0.01 


Core 






20 ±4.0 


0.55 ±0.11 


14±11 


0.64±0.13 


C7 






14±2.1 


2.82 ±0.24 


14±5 


1.68 ±0.34 


C3 


1994.70 


5.0 


284 ±42.6 






0.18 ±0.04 


Core 






26 ±3.9 


■ 


• 


0.21 ±0.04 


C7 






13 ± 1.9 


1.86±0.17 


18±5 


0.77 ±0.15 


C5 






12± 1.8 


3.00 ±0.17 


15 ±3 


0.72 ±0.14 


C3 






7± 1.0 


4.34 ±0.32 


8±4 


1.45 ±0.29 


C2 


1995.15 


22.2 


724 ±108.6 


_ 


_ 


0.07 ±0.01 


Core 






24 ±3.6 


0.23 ±0.04 


12±9 


0.07 ±0.01 


C8 






9 ±1.3 


1.17 ±0.06 


12±3 


0.07 ±0.01 


C6 






2 ±0.3 


2.37 ±0.19 


14±4 


0.15 ±0.03 


C4 






10 ±1.5 


3.43 ±0.10 


15±1 


0.20 ±0.04 


C3 


1995.31 


22.2 


421 ±63.1 


_ 


_ 


0.10±0.02 


Core 






16 ±2.4 


0.27 ±0.06 


12± 11 


0.10 ±0.02 


C8 






7± 1.0 


0.77 ±0.08 


14±6 


0.10±0.02 


C7 


1995.47 


22.2 


174 ±26.1 






0.04 ±0.01 


Core 






14 ±2.1 


0.31 ±0.05 


29±9 


0.07 ±0.01 


C8 






6 ±0.9 


0.85 ±0.07 


18±4 


0.06 ±0.01 


C7 






6 ±0.9 


2.07 ±0.19 


23±5 


0.43 ±0.09 


C5 


1995.65 


22.2 


220 ±33.0 






0.07 ±0.01 


Core 






81 ±12.1 


0.16 ±0.02 


10±6 


0.05 ±0.01 


C9 






29 ±4.3 


0.40 ±0.07 


13 ±10 


0.32 ±0.06 


C8 


1995.65 


8.4 


252 ±37.8 






0.11 ±0.02 


Core 






32 ±4.8 


0.43 ±0.10 


9 ±12 


0.10 ±0.02 


C8 






8± 1.2 


1.36±0.18 


15±7 


0.53 ±0.11 


C6 






3 ±0.4 


2.64 ±0.26 


20±5 


0.40 ±0.08 


C4 






5 ±0.8 


3.97 ±0.34 


15±5 


1.19±0.24 


C3 


1996.34 


22.2 


221 ±33.1 


_ 


_ 


0.08 ±0.02 


Core 






35 ±5.3 


0.43 ±0.06 


9±7 


0.24 ±0.05 


C9 






10 ±1.5 


0.97 ±0.09 


12±5 


0.15 ±0.03 


C7 






6 ±0.9 


1.47±0.13 


5±5 


0.20 ±0.04 


C6 


1996.53 


15.3 


225 ±33.8 


_ 


_ 


0.08 ±0.02 


Core 






17 ±2.6 


0.33 ±0.05 


9±8 


0.08 ±0.02 


C9 






15 ±2.2 


0.75 ±0.09 


4±6 


0.23 ±0.05 


C8 






4 ±0.6 


1.33 ±0.23 


17±9 


0.43 ±0.09 


C7 






3 ±0.4 


4. 12 ±0.34 


15±4 


0.71 ±0.14 


C3 






1±0.1 


5.14 ±0.24 


27±2 


0.12 ±0.02 


C2 


1996.60 


22.2 


218 ±32.7 


_ 


_ 


0.06 ±0.01 


Core 






17 ± 2.6 


0.73 ±0.09 


6±6 


0.26 ±0.05 


C8 






3 ±0.4 


1.16±0.19 


15±9 


0.22 ±0.04 


C7 






4 ±0.6 


1.78 ±0.08 


19±2 


0.05 ±0.01 


C6 


1996.63 


5.0 


183 ±27.4 






0.03 ±0.01 


Core 






22 ±3.3 


1.19 ±0.07 


9±3 


0.24 ±0.05 


C7 






8 ±1.2 


4.09 ±0.27 


20±3 


1.19±0.24 


C3 






7 ±1.0 


9.58 ±0.49 


23±2 


3.39 ±0.68 


CI 






3 ±0.4 


12.06 ±0.44 


10±2 


1.14±0.23 


X 


1996.82 


15.3 


251 +37.6 






0.01+0.00 


Core 






13 + 1.9 


0.40 ±0.06 


19±7 


0.08 + 0.02 


C9 






10± 1.5 


1.49±0.16 


12±6 


0.51 +0.10 


C7 


1996.90 


22.2 


260 ±39.0 






0.04 + 0.01 


Core 






24 ±3.6 


0.22 ±0.04 


22± 10 


0.08 + 0.02 


CIO 






3 ±0.4 


0.76 ±0.12 


16 ±8 


0.08 + 0.02 


C8 






7 ±1.0 


1.54±0.16 


7±6 


0.38 + 0.08 


C7 






3 ±0.4 


4.69 ±0.29 


13±3 


0.50 ±0.10 


C3 
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Table 2. Modelfit parameters (continued from previous Col.). 



Epoch 


V 

[GHz] 


r mjy . 
L beam ^ 


r 

[mas] 


'P 
[°] 


e 

[mas] 


Id." 


1997.03 


8.4 


194 ±29.2 


- 


— 


0.09 ± 0.02 


Core 






22 ±3.3 


0.77 ±0.07 


12 ±5 


0.23 ± 0.05 


C8 


1997.58 


22.2 


929 ± 139.4 






0.07 ±0.01 


Core 






6 ±0.9 


0.24 ±0.14 


■ ■ 


0.01 ±0.00 


CIO 






32 ±4.8 


0.51 ±0.07 


7±7 


0.32 ± 0.06 


C9 






7± 1.0 


1.69 ±0.07 


11 ±2 


0.06 ±0.01 


C7 


1997.93 


8.4 


385 ± 57.7 






0.10 ±0.02 


Core 






19±2.8 




■ 


0.16 ±0.03 


C9 






15 ±2.2 


■ 


7 ± 3 


0.30 ± 0.06 


C7 






7 ± 1.0 


3.18 ±0.24 


15 ±4 


0.84±0.17 


C5 


1999.41 


8.4 


922 ± 138.3 


- 


- 


0.10 ±0.02 


Core 






60 ± 9.0 


0.68 ± 0.07 


10±5 


0.10 + 0.02 


CIO 






17 ± 2.6 


1.58 + 0.08 


13 + 2 


0.20 + 0.04 


C8 






4 ±0.6 


3.16±0.19 


6±3 


0.28 ± 0.06 


C6 






4 ±0.6 


3.78 ±0.22 


15 ±3 


0.40 ±0.08 


C5 


1999.55 


15.3 


1 157 ± 173.6 






0.06 ±0.01 


Core 






54 ±8.1 


0.40 ± 0.03 


15 ± 3 


0.08 ±0.02 


CI 1 






26 ± 3.9 


1.12±0.09 


12±4 


0.39 ±0.08 


C9 






2 ±0.3 


1.89 ±0.26 


10 ±8 


0.28 ± 0.06 


C8 






10 ± 1.5 


2.90 ± 0.24 


9±4 


1.17±0.23 


C6 


1999.89 


5.0 


558 ±83.7 


- 


- 


0.18 ±0.04 


Core 






39 ±5.9 


1.29±0.11 


8±5 


0.13 ±0.03 


C9 






19 ±2.8 


2.64 ±0.15 


10 ±3 


0.83 ±0.17 


C7 






6 ±0.9 


4.43 ±0.29 


12 ±3 


1.00 ±0.20 


C4 






1 ±0.1 


6.92 ±0.98 


18±8 


1.93 ±0.39 


C2 






7± 1.0 


11.86 ±0.58 


17±2 


4.79 ±0.96 


CI 


2000.82 


5.0 


520 ±78.1 


- 


- 


0.12±0.02 


Core 






35 ±5.3 


0.77 ±0.12 


20 ±8 


0.10 ±0.02 


Cll 








1 A'i 4- n no 


12 ± 3 


U.^Z ± U.Uo 


C9 






11 ±1.6 


2.71 ±0.18 


11±3 


0.74 ±0.15 


C7 






9 ±1.4 


4. 16 ±0.26 


11±3 


1.25 ±0.25 


C5 






8 ±1.3 


11.16±0.51 


17±2 


4.43 ±0.89 


CI 


2001.17 


15.3 


549 ±82.4 






0.06 ±0.01 


Core 






72 ±10.8 


0.45 ±0.04 


20±4 


0.18±0.04 


C12 






14±2.1 


1.81 ±0.16 


17±5 


0.70 ±0.14 


C9 






5 ±0.8 


3.06 ±0.22 


13±4 


0.49 ±0.10 


C7 






3 ±0.4 


5.08 ±0.29 


16±3 


0.51 ±0.10 


C4 



" Identification of the individual components. If a component ap- 
peared only in a single epoch it is labelled with X. 



